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We consider the Casimir effect of a massive vector fieid between two semi-infinite dielectric slabs. 
We first derive the generalization of the Lifshitz formula that gives the Casimir interaction energy 
of two magnetodielectric slabs separated by a magnetodielectric medium due to the vacuum fluctu- 
ations of a massive vector field. We then discuss the asymptotic behaviors of the Casimir energy 
and the Casimir force in various limits, such as low temperature, high temperature, small mass, 
large mass, up to the first order in the finite conductivity correction, for two real metal semispaces 
s ! ■ whose dielectric property is described by the plasma model. Application to the Casimir effect in 

I Randall- Sundrum spacetime is briefly discussed. 

<N 

. I. INTRODUCTION 

o 

C^l ■ The Casimir effect is an interesting phenomena that is due to the vacuum fluctuations of quantum fields. The 
Casimir effect of massless scalar field, electromagnetic field (masslcss vector fields) and massive scalar field have been 
extensively studied (see e.g. [IH1]). In contrast, massive vector field has not been much considered, mostly due to 
O J its substantial complication compared to the massless case. Nevertheless, massive vector field plays an important 
role in the study of physics with extra dimensions. To study a quantum field in a spacetime with extra dimensions, 
c~j ■ one can use the Kaluza-Klein decomposition to decompose the quantum field to an infinite tower of massive fields 
in four dimensions. This approach has been intrinsically used in a number of works to study the Casimir effect of 
scalar field or spinor field in spacetime with extra dimensions such as Kaluza-Klein spacetime and Randall-Sundrum 
^ , ' spacetime. For electromagnetic field, thanks to the work of Barton and Dombey the Casimir effect between two 
parallel perfectly conducting plates in Kaluza-Klein spacetime of the form M x S 1 and in Randall-Sundrum model 
have been computed respectively in Q and Q. 

The Casimir effect of a massive vector field was first considered in and an explicit formula for the Casimir 
energy has been derived for two perfectly conducting plates of finite thicknesses. Asymptotic behaviors of the massive 
corrections were computed. Unlike the electromagnetic field, a massive vector field has three polarizations. Besides 
two transverse polarizations that correspond to the two polarizations of an electromagnetic field, a massive vector 
field also have longitudinal modes and these latter modes can penetrate through perfectly conducting objects. In a 
recent work Q, we generalized the work Q to magnetodielectric slabs. We found that for Casimir effect of massive 
vector field on general magnetodielectric slabs, one cannot separate the contribution from the transverse magnetic 
modes from the contribution of longitudinal modes, and we wound up with a very complicated formula for the Casimir 
interaction energy. The TM contribution, which is a combination of the contributions from the transverse modes and 
the longitudinal modes, is expressed in terms of log determinant of a four by four matrix. It is almost impossible to 
■ study such a formula analytically. 
$H ' To study the electromagnetic Casimir effect of two semi-infinite dielectric slabs, Lifshitz Q has derived a formula 
which was later named after him. The first goal of this article is to derive a formula that looks more similar to the 
Lifshitz formula for the Casimir interaction between two semi-infinite magnetodielectric slabs due to the fluctuation 
of a massive vector field. The second objective is to consider the deviation from a perfect conductor. As in the works 
[lol fill ] , we model the two semi- infinite slabs as real metals whose dielectric property is described by plasma model: 

e(w) 

where £0 is the vacuum permittivity and oj p is the effective plasma frequency. We expand the generalized Lifshitz 
formula perturbatively in the dimensionless variable 

c 
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where d is the distance between the slabs, a — gives the perfect conductor limit, and a small a amounts to finite 
conductivity corrections. For the zero and first order in a, we study the asymptotic behaviors of the Casimir energy 
and the Casimir force in various regions: large mass high temperature, large mass low temperature, small mass high 
temperature and small mass low temperature. In the last section, we also briefly discuss the application of these 
results to the Casimir effect of real metals in Randall-Sundrum spacetime. 



II. MASSIVE VECTOR FIELD 



In this section, we recall some basic facts about massive vector fields that are required for the study of Casimir 
effect. For more details, one can refer to [H, H|. 



A. Proca equations 



Define the electric field E and the magnetic field B by 

BA 

E = — r--V<£, B = VxA. (1) 

dt w 

The Proca equations for a massive vector field propagating in a medium with permittivity e and permeability \x are 

V B = 0, 

VxE + 5 = 0, 
at 

m 2 ± (2) 
V-D + — <S> = Ph 

„ TT 9D m 2 c 2 . 

VxH -^ + 7^ A = J/ - 

The first two are the well-known Maxwell's equations which are automatically satisfied because of ([1]). The continuity 
equation 

t + V- a ,=0 (3) 



implies that the Lorentz condition 



has to be satisfied. 
Let 



and assume the linear relations 



/oo 
duj<fi(x., to) 
-oo 

/oo 
dwA(x,w)e-^*, 
-oo 

D(x,w) = e(w)E(x,o;), 
H(x,w) = -^tB(x,u). 

When there are no free charges (i.e. p/ = 0) and free current (i.e. J/ =0), the Proca equations ([2]) are equivalent to 
the following two equations for 0(x, u) and A(x, u): 



„ uj" in 



V x V x A - e/ic 2 V(V • A) - ( e/iw 2 — ) A =0. 



(5) 
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The Lorcntz condition ((4]) becomes 

<£=-— V-A. (6) 
The solutions of these equations can be divided into two types: 



• 2 

IC 



• Transverse waves. These are waves with V ■ A = 0, <f> = and 

V x V x A - ( efiuj 2 - ^-f- ] A = 



2 2 

m c 

To solve this, we look for a vector wave function A satisfying 

VxVxA=4 

Then lo satisfies the dispersion relation 



7 2 2 
%. = £[LLO — 



h 2 



The solutions can be divided into two families A TE and A™, which in the massless limit correspond to the TE 
(transverse electric) and TM (transverse magnetic) waves of electromagnetic fields. 

Longitudinal waves. These are waves with V x A = 0. The equations ([5]) become 

2 2 \ 



efic 2 V(V ■ A) - ( efico 2 p- ) A = 0. 



To solve this system, one first find a scalar function cf> satisfying 
Then A is given by 



-V 2 = k 2 L cj>. 



V0, 



k 2 c 2 



and the dispersion relation is 



In this case, B = and 



9 9 

L c 2 e[iY? 



H 2 efik 2 

B. Plane waves solutions 

In this article, we only consider plane waves. Choosing z as the distinctive direction, the plane waves are 
parametrized by kj^ = (k x , k y ) £ M. 2 . The transverse TE and TM waves are: 



TF reg 1 / r~o 7" 

A ki ' out (x,w) = — V X ^ e ^x+ l k v y T t y /k- T -k 1 z ez 

k_ 

TM,5, % 1 



e^ e^, | e 

k± k± 



A 



' out 



k. 



(x,Cj) = V X V X ( e ik**+ik V VTW k T- k i z e z ) 

kxk± V / 

fca; \/ fc^ k 2 k v \/ k^p k\ kj_ 

± — )L r L 1 ^-e, ± av 1 ±e y + — ( 

ferKi krk± kx 



ikxZ+ik y y=fiy/k%. — k^_z 
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The longitudinal waves are given by 



j rog r 2 fr 2 

ut (x Ul) — L gikgX+ikyy^i^ k\ ~k 2 , z 



'k 



Here reg and out stand for regular waves and outgoing waves respectively. 



III. FINITE TEMPERATURE CASIMIR INTERACTION OF TWO MAGNETODIELECTRIC SLABS 

SEPARATED BY A MAGNETODIELECTRIC MEDIUM 




* 

FIG. 1: Two parallel magnetodielectric slabs separated by a magnetodielectric medium. 



Consider two parallel slabs, each of area A, occupying the regions — L < z < and d < z < L, with permittivities 
£z,e r and permeabilities fii,fi r respectively. The medium between the slabs is assumed to have permittivity e& and 
permeability (see Fig. [I}. As discussed in [j| the boundary conditions of a massive vector field are the 
continuities of A, /i _1 B|| and the continuity of either (f> or <9„A n , where n is a unit normal vector to the boundary. 

In the following, we compute the contribution to the Casimir free interaction energy from the transverse TE modes, 
and from the combination of transverse TM modes and longitudinal modes. We denote the former as TE contribution, 
and the latter as TM contribution. 



A. TE contribution 



For TE contribution to the Casimir free interaction energy, </>(x, t) = and 

A M =ApJu J°°J^ ^^(a,(kx)A^/«(x, W )+«8.(kx)A™ : r(x,«)) ^ ■ 
Here * = I when —L < z < 0, * = b when < z < d and * = r when d < z < L, 

Pt.* = J k\ ~ — k\. Taking the limit L — >• oo amounts to setting *B; = and 2l r = 0. The continuities of A and 
/U _1 B|| then give 

21; = a 6 + Q5 6 , 
Mi Mfc 

_ PliLfB re iPT,rd = Pr^ (% be - tP T,„d _ ^ be WT, b d\ _ 

Mr Mb 
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Eliminating 21/ and *8 r from these equations, we obtain 

(PT,bPl +PTJ^b) ( Bb = (PT.bfM - PT,lHb)%b, 

(p T ,bPr ~ PT^b) < B b e ip ^ d = (p T ,b»r+PT,r»b) Slfte"^*"*. 

To have nontrivial solutions for (Sit,, 23&), we find that the eigenfrequencies ui should satisfy the equation 

(PTMPl +PT,tfJ-b) (PT,bPr + PT,r^b) e~ WT - bd - (pT,bM ~ PT,l^b) {PT,bVr ~ PT.rPb) e WT - bd = 0. 

Factoring out 

(PT,bHl + PT,lPb) (PT,b^r +PT,rPb) e~ WT - bd 

which can be attributed to the self-energies of the individual slabs, we obtain 

1 _ (PTMM - PT,lVb) iPT,bHr ~ PT,rPb) e 2ip T , b d = q (g) 
{PT,b^l +PT,lHb) (PT,bHr +PT,rHb) 

Alternatively, we can derive this by writing ([7]) in the form 



where 

R TE = 

PT,*/J-b + PT,b[i. 



D TE PT.*Hb - PT,bH* , 

R t = ■ , * = l or r. 



Then 

93 6 = -Rj E K b = Rj E Rj E e 2 ^< bd <& b . 

From this, we find that there are nontrivial solutions for 55/, if and only if uj satisfies the relation ([5]). Hence, by using 
the standard contour integration technique and the Matsubara formalism, we find that the TE contribution to the 
Casimir free interaction energy is 



^Cas — 



k R TA c 



2tt 
where 

2T7nksT 

are the Matsubara frequencies, 



V' / dk ± k ± In (l - Rr^n)Rj E ^n)e- 2qTA ^ )d ) , (9) 



R * ( l = m cc\ , 7F\ — F*v * = 1 or r ' 10 ) 



q T ,*(0 = \J s,(iO^*(^)C 2 + -jp- + * = l,bovr. 

B. TM contribution 

For the TM contribution to the Casimir energy, which comes from the transverse TM modes and the longitudinal 
modes, let 
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and 

-*£*££££ 

x (e,(kx)A™;'«(x,u,) + £,(k ± )A™;"(x.") + e.fkxJA^fx.a.) + S,(k i )A^'(x, w ) 
where * = I when — L < z < 0, * = b when < z < d and * = r when d < z < L: 

T rc S C^fc? 
A 'out(„ , ,\ _ L >* ik x x+ik v y^ip L .,z 

^• A ' >~k ± LO e 



A^°f (x, W ) = (ge* + ^e, T f^e z ) e **.*+*.«fF*i,.*. 



Here 



PL,* = Jki - kj_. 



Taking the limit L — > oo amounts to setting 55/ = 0, Si = 0, £ r = and <£ r = 0. The continuities of A, and /i 1 B|| 
then give the following two sets of relations: 

pic, + gj = Pl± {€b _ Sfc) + (g 6 + 

^ - ^ = + Jfc) - ^ - ?6), 

- ^>2) r e ipT '" d + ff r e* x " rd = ^-(€ b e- ipT - bd - ® b e ipT < bd ) + (C 6 e _ * t '» d + &,e <pi -<* d ), 
Eliminating (£;,2) r ,3V give relations of the form 

®b\ _ nTM ( 

o e -iPL, b dj y£ b j n r y o gipi.i.dj ^ b 

where i2 z and i? r are 2x2 matrices. Hence, the TM contribution to the Casimir free interaction energy is 



TE _ k B TA 



V' / dk ± k± lndet (1 - i2™(ien)C/(^n)«™«r I )C/(^n) ) , 



where 



2tt 

n=0 



' pTM nTM ' 

= r> , -K* = nTM nTM , * — i Ol T. 
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For * = I or r, 



R TM 



nTM 



=1 1 — : ! 



2 «T *lL,bKT,bVb ™ 2 " 2 - 2 



L.b 

m"c" n 



A* ft 2 



"6"* 



nTM _ 2 lT,bK 2 L ^k 2 ± 2 2 
■ft*, 21 — °* a { K T.b^* ~ K T.*^b), 



TM 
22 



2«i *9i,6 [qT,*fJ-*KT b + 9T,6M6«T* 
1 ■ : 



A„ 



1 for * = I and cr* = — 1 for * = r. For * = I, b, r, 



k t ,* =\j e*(i£)/i*(i£)£ 2 H — ^2—, <?t,* = y K T,* + fc l' 



c 2 ' h 2 e*{i£,)n*{i£,)' 



This presentation of the TM contribution has substantially simplifies the presentation given in j8[. Moreover, it has 
been cast in the form of TGTG formula [HI, [l4| , which is the general presentation for the Casimir interaction energy 
between two objects. 



IV. MASSLESS LIMIT 



In this section, we consider the limiting case when the mass is zero. When m = 0, it follows immediately from (|9]) 
that the TE contribution to the Casimir free interaction energy is given by the same formula (j9j), with Rj E given by 
(TTU|) . but q Tt: is reduced to 



<?7>(£) = ye^OM^X 2 + fc L * = Z,6orr. 
For the TM contribution, obviously, when m = 0, 



nTM _ nTM _ n 
^.12 — -fVl2 — u ' 



On the other hand, when m = 0, 



Therefore, 



kh.i = k L j- = k L , b = - q Lt i = q L . r = q L , b = d — + k\ 



nTM _ nTM _ n 
jxj 99 — n„ 99 — u. 



4,22 — • n r,22 

Hence, the TM contribution to the Casimir free interaction energy reduce to 



n=0 

where 



=^ £' f dk^k ± in (1 - ^( 4 e„)i?™(^j e - 2 ^ (? " )d 



qT,*^*£b^b + qT,b^b£*^* 

qr,*£b — qr,b£* 



<7t,*£& + qr,b£* 

These recovers the Lifshitz's formula [i| for the electromagnetic Casimir free interaction energy between two dielectric 
slabs. 
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V. THE LIMITS OF PARALLEL PERFECT CONDUCTORS SEPARATED BY VACUUM 



In this section, we consider the case where two perfectly conducting plates are separated by vacuum. This can be 
achieved by taking the limit Si,s r — ¥ oo and setting e b = e , /i^ = Ho- I n this case, we find that 



«T,fc =KL,b = 



qr,b =qL,b 





m 2 c 2 


c 2 






9 9 




h 2 



Kb, 



k\ := q b , 



for * = / or r, 



c 



pTM 

pTM 
H *,12 



pTM , _ 



£ 2 
C z 



3o 

2 i „ ,.2 



i 2 c 2 fc 2 



-qb(qon b + gb^p) h — p-^ _ -<?o(<?6 + fci) + 2g fe fc 2 

9b ("Zo^ 2 + % Kq) 



2gb'«b I %r 



m 2 c 2 fc 2 



9b(9o« 2 H p- 



9o (3? - fcl) 
2q b K b (q b - g ) 



h'b 



9h(<?oK 2 + 9b«o) 

2 „ 2 i, 2 



m 2 e 2 fc 2 



fi 2 



2g b (g 2 - fc 2 )fcj 
" K b q {qb + qo)(q b -kj)'' 



R 



TM 



9b(9o« 2 - q b n 2 ) H p-^ _qb-qo 9o(9 2 + fcj.) + 2g b fcl 



/ 2 i 2 \ i "i 2 e 2 fc 2 i 



qb + qo qoiql-li]) 



(11) 



One can then show that 



(12) 



and 



det (l - R™UR™U) -> (1 - e~ 2 ^) (l - A 2 e " 2 ^) , 



(13) 



where 



A = 



qb - go 
qb + qo' 



For the zero Matsubara frequency io, we have to use the Schwinger-DeRaad-Milton prescription [15( where the limit 
£l,e r — > oo is taken before setting £o = 0. 

From (fT2")l and ([TB")) . we find that in the perfect conductor limit, the Casimir free interaction energy is given by 



; =-7; X/' / *-L^J- ln 



2tt 



-n JO 



ri=0 



1 - e 



-2<ft(f„)<i 



9b(£n) - 9o(£n) 



qb {in) + 9o(£n) 



3 -2 9t (£„)d 



In the language of 0], the factor (l - e ^ 2qb ^ n)d 



gives two discrete mode contributions, whereas the factor 



1 - 



qb {in) - qo{in) 



-1 2 



96 {in) + qo {in) 



3 -2 96 ({„)d 



gives a continuous mode contribution. 
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VI. ASYMPTOTIC BEHAVIOR FOR REAL METALS DESCRIBED BY PLASMA MODEL 

In this section, we consider the case where the two semi-infinite slabs are metals with permittivities described by 
plasma model: 

and permeabilities fj,i = /x r = (j®. Here lo p is the effective plasma frequency. 

The medium between the metals is assumed to be vacuum, i.e., Eb = So, /U& = Mo- Let 

mcd c a p £ n d 

A=_ ^~' a= — 3 = 7T~7: z n = • 

n uivd Zira c 



Then 



so 

Making a change of variables 



we find that the TE Casimir free interaction energy is given by 

^as=^£' /^dzzln(l- [i? TE ]%~-), (14) 

where 



n=0 • / V Z "+ A 



TE VI + a 2 z 2 - az 
-K = 



vT-j- a 2 z 2 + az 
The TM Casimir free interaction energy is more complicated: 



l=0 W^+ x 

where 



« =^ £' /T <^ det f 1 - We" 22 ) , (15) 



TM 1 / Vil G*V\2 



* U \-<r*V 21 V 22 

U = (a^l + z 2 a 2 [zl + A 2 ] + z [l + a 2 (z 2 + A 2 )]) (zz 2 [l + a 2 (z 2 + A 2 )] + + a 2 z%y/z 2 - A 2 + a 2 z 2 z 2 [z 2 + A 2 ]) 
+ A 2 (z 2 -z 2 -A 2 ), 

= (a^l + ^a 2 ^ + A 2 ] - z [l + a 2 (z 2 + A 2 )]) (zz 2 [l + a 2 (z 2 + A 2 )] + \/l + a 2 z 2 ^z 2 - A 2 + a 2 z 2 z 2 [z 2 + A 2 ]) 
+ A 2 (z 2 -z 2 -A 2 ), 



V 12 =2A 2 zv/z 2 + A 2 (1 + a 2 [z 2 + A 2 ]) , 
_ 2zz 2 (lW[z 2 +A 2 ])(z 2 -z 2 -A 2 ) 

V^T+a 2 

y 22 = (aVl + ^a 2 ^ 2 + A 2 ] + z [1 + a 2 (z 2 + A 2 )]) (-zz 2 [l + a 2 (z 2 + A 2 )] + y/l + o?z 2 Jz 2 - A 2 + a 2 z 2 z 2 [z 2 + A 2 ] 
+ A 2 (z 2 -z 2 -A 2 ). 
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A straightforward computation gives 

dot (l - Rmre^) =1 - m±M±^i e -- + iVuVv + VM 2 

= (l-T +e - 2z ) (l-T_e" 2zN ) 



e 4z 



where 



Hence, 



_ VA + V22 + 2^12^21 ± (V22 - Vnh/^i + V22) 2 + 4Vj 2 y~ 
± 2U 2 



771TM _ksTA 

^Cas 



27Td 2 

n=0 



/ dzzln[(l-T + e- 2 *)(l-T.e- 2 *)]. (16) 



As in [ijj, [TT( j we consider the asymptotic behaviors of the Casimir free interaction energy when a -C 1 ■ a = 
corresponds to the perfect conductor limit. Nonzero a corresponds to finite conductivity corrections. For aluminium, 
gold and copper plates, the plasma wavelength a p are given respectively by 98nm, 132nm and 132nm [l6[. With 
separation between the plates d in the range lOOnm = 0.1/xm to 1/xm, a < 0.25. 

The computations of the asymptotic behaviors are quite tedious and we leave them to the Appendix [A] Here we 
discuss the results. 

Let 

27tkgTd ttT 

He T e g ' 

where the effective temperature 

he 



T, 



eS 



2k B d 



is defined in [17]. 

In the large mass high temperature A, $ ^> 1 region, we find that 

3k B TA^ _ 2X ( ^ 8A 



^~-^Ae-"ll- T a + ...]. (17) 



Differentiating with respect to d, we find that the Casimir free interaction force 

p _ dEcas 
-^Cas ^ , 

ad 

behaves like 



^Cas ^V 2 Ml--«+... . (18) 



And 3 ' V 3 

Here for each order in a, we only listed the leading term in A. Observe that these terms are linear in temperature, 
coming from zero Matsubara frequency, but exponentially small in A. 
In the large mass low temperature -d <C 1 <C A region, 



3hcA , 3 2A ( , 8A 
ZhcA , 5 , A / 8A 



(19) 



Similarly, for each order in a, we only listed the leading term in A. This terms comes from the zero temperature 
region. They are exponentially small in A. The thermal corrections are exponentially small. 

From (|18p and (|19[) . we find that in the large mass region, a small deviation from infinite conductivity tends to 
reduce the strength of the Casimir force. 
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Next, we turn to the small mass A <C 1 region. In the high temperature ■& S> 1 region, 



E Ca 



F, 



Cas 



k B TA 
8nd 2 



k B TA 
And 3 



Cfl(3) 



Cfl(3) 



1 

2 
1 - 



Cfl(3) 



1 



A 2 In A 



-A 2 + 



(21n2-l)A 2 



2Cfl(3) 
4 



Cfl(3) 



A 3 + 



Cfl(3) 3Ci?(3) 



A' 5 + 



(20) 



1 



-A 2 + 



1 



Cfl(3) C«(3) 



A d + 



-3 



1 - 



3Cfl(3)' 



Again, these terms are linear in temperature, coming from zero Matsubara frequency. In the masslcss A = limit, we 
find that 



Eca 



k B TA 

8ird 2 
k B TA 

And 3 



C«(3){1- 2a +...}, 
C«(3){1- 3a + ...}. 



(21) 



Setting a = gives 



-Etas 

Fcas 



k B TA 

8ird 2 
k B TA 

And 3 



Cr(3), 
Cfl(3), 



(22) 



which are well known results for high temperature asymptotics for perfect conductors 0, [H, [l9| . From (|2"0]l , it is 
interesting to observe that the existence of small mass and small deviation from infinite conductivity both reduce the 
strength of the Casimir force. Moreover, the existence of small mass will reduce the effect of small deviation from 
infinite conductivity. 

In the small mass low temperature A, $ <C 1 region, we consider two cases. When < A < 1, we find that 



Ecas — — 

Fcas 



it 2 he A 
720d 3 

TT 2 hcA 

240d 4 



1 - -^A 2 + 
1-^A 2 + 



/90 


5) 






a + . . .1 




-4 


[i-V 










7T^ 



16 








i - -V" 


a+ ...| 


~ y 


2tt 2 





These terms come from the zero temperature region. The thermal correction terms are exponentially small. 
When A <C i? <C 1, there are polynomial order terms of the thermal correction: 



Cas 



it 2 he A 
720d 3 



1 



15 



A 1 



60 



lni?-lnA-ln(27r) + 



-4 



1 



2- 



,-A 2 1 



2- 



¥ Cfl(3)^ 



Cas 



ir 2 hcA 
240d 4 



1 tf 4 

+ 3^ 

16 " 

~y 

In the masslcss limit A = 0, we have 

ir 2 hcA 

Ecas = — 
Fcas = — 



1-SCh(3)^ 3 ~^A 2 



8tt 6 



4tt 4 



720d 3 
n 2 hcA 
240d 4 



^C,(3)^ 3 - J ) 



7T 

ltf 4 
3^ 



45 
2^6 



C«(3)^ 3 



^4 



16 

T 



45 
8^ 



Cfl(3)^ d a 



(23) 



The latter agrees with jllj up to the first order in a. Again, we find that the existence of small mass and small 
deviation from infinite conductivity both reduce the strength of the Casimir force. Moreover, the existence of small 
mass will reduce the effect of small deviation from infinite conductivity. 
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VII. APPLICATION TO CASIMIR EFFECT IN RANDALL-SUNDRUM SPACETIME 



One of the main motivation for studying Casimir effect of massive vector field comes from extra-dimensional 
physics. To study the Casimir effect of the electromagnetic field in a spacctimc with extra dimensions, one can use 
the Kaluza-Klcin decomposition to decompose the electromagnetic field to an infinite tower of massive vector fields 
in four dimensions. For two parallel perfectly conducting plates in the Randall-Sundrum spacctimc, we have studied 
the Casimir effect along this line in [7j. 

Recall that the spacetime underlying the Randall-Sundrum (RS) model is a hD anti-de Sitter space (AdSs) with 
background metric 

ds 2 = e- 2K ^i lab dx a dx b - dy 2 , (24) 

where rj a b = diag(l, —1,-1,-1) is the usual AD metric on the Minkowski spacetime M 4 . The extra dimension with 
coordinate y is compactified on the orbifold S 1 /Z2. K is a parameter that determines the degree of curvature of the 
AdSs space. There are two 3-branes with equal and opposite tensions, localized at y = and y — -kR respectively, 
where R is the compactification radius of the extra dimension. Using Kaluza-Klcin decomposition, an electromagnetic 
field in the 5D Randall-Sundrum spacetime is decomposed into an electromagnetic field (massless vector field) in the 
AD Minkowski spacetime, and an infinite family of massive Proca fields with masses mi < mi < . . ., where 



rrijC 



are positive solutions of the equation @, H3] : 



/ ZG^ \ / Z \ f Z&^^^ 



Jo {~J Yo {—)- Yo [-J 

To solve the hierarchy problem between the Planck and electroweak scales, it is required that kR ~ 12 (20|. In this 
case, one can show that 

n 

For the effect of the extra dimension to be significant, we require that 

irnde- 127 ' ~ 1. 

For d ~ lOOnm, this amounts to k ~ 10 8 GeV, consistent with the numerical result found in Q for perfect conductors. 
For k > 10 10 GeV, 

rriicd 
A; = ^-»l, 

we find from the asymptotic analysis in the last section that the leading term of the Casimir force due to the Kaluza- 
Klein nonzero modes is exponentially small. Hence we only see the AD (massless) electromagnetic Casimir effect. 
When k < 10 8 GeV, there will be significant correction to the Casimir force from a finite number of the Kaluza-Klein 
nonzero modes with 

mjcd 
A ' = -h <L 

To get a rough idea of how big the contribution can be, one can take a finite sum over the asymptotic formulas derived 
in the previous section. 



VIII. CONCLUSION 



In this article, we have derived a generalization of the Lifshitz formula to describe the Casimir interaction between 
two parallel semi-infinite magnetodielectric slabs separated by a magnetodielectric medium. Specialized to two parallel 
real metals with dielectric property described by the plasma model, we expand the Casimir energy and the Casimir 
force up to first order in the finite conductivity correction. Asymptotic behaviors of the Casimir energy and Casimir 
force in combinations of high/low temperature and large/small mass regions are discussed. It is found that when the 
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mass is large, the Casimir force is exponentially small. In the small mass region, the dominating term is the masslcss 
term, and the massive correction tend to reduce the strength of the Casimir force and the effect of deviation from 
infinite conductivity. The application to Casimir effect in Randall-Sundrum spacetime is briefly discussed. It is found 
that for the electromagnetic Casimir effect in the Randall-Sundrum spacetime to be significantly different from the 
4D electromagnetic Casimir effect, the warping parameter of the Randall-Sundrum spacetime n must satisfy n ~ 10 8 
GeV, when the separation between the two metals is about lOOnm. 

A massive vector field has longitudinal modes which satisfy dispersion relation different from the dispersion relation 
of the transverse modes. However, the contribution of the transverse magnetic modes and the longitudes modes 
cannot be separated. This phenomena has also been observed in the case of two concentric spherical bodies [2l| . 
This makes the Casimir effect of a massive vector field considerably more complicated than the Casimir effect of an 
electromagnetic field. In this work, we have only expanded the Casimir energy up to the first order in the finite 
conductivity correction. However, our scheme can be used to obtain higher order conductivity corrections. 
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Appendix A: Asymptotic behaviors of the Casimir energy 



In this section, we compute the asymptotic behaviors of the Casimir free interaction energy in various limits. We 
consider the TE and TM contribution separately. 



1. TE contribution 



Expanding the logarithm in (1141) . we find that 

oo oo 1 „oo 



Since 



when a <C 1, we have 



where 



7 rp a ^ ^ poo 

£TE = _ kBlA y,yl dzz F^TE] * e - 2jz 



[i? TE ] 2j =l-4ajz+..., 

1 T 1 A 00 °° 1 poo 

WA^,yl / dzze -2j* 



n=0 j=l- 



A straightforward integration gives 



(Al) 



71=0 j = l 



2nd 2 



' 2(4 + a 2 ) 2^iTy A . 2jV ^ 



(A3) 
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Let 



$ = 



He 



so that z n = n§. $ 3> 1 corresponds to high temperature limit, and d <C 1 corresponds to low temperature limit. 

In the large mass A > 1 limits, obviously Ao and Ai goes to zero exponentially fast. In addition, if $ ^> 1, the 
leading terms come from the terms with zero Matsubara frequency. From (| A2|) and (|A3|) . we find that in the limit 

k B TA 2 _ 2A 

When i9 <C 1 <C A, the dominating term comes from the zero temperature term, and the thermal correction terms are 
exponentially small. We have 



J 

hcA 



Alt 2 d? 

HcA 

47T 2 d 3 
47T 2 d 3 



dzzy/ z 2 — X 2 e 



dz(z + \)y / z 2 + 2zXe-' 2{z+x) 

V2\h~ 2X r 'dzzh- 2z 

Jo 



-A3 



,-2A 



167T 2 (i 3 

Similarly, we find that when ^ < H A, 

. HcA 5 o \ 
A § — A5e" 2A . 

47T2d 3 

Next we consider the small mass A« 1 limits. In the high temperature d 3> 1 region, the terms with nonzero 
Matsubara frequencies go to zero exponentially fast. Using the formula 



_ z _ 1 

we find that the asymptotic behavior of Ao is given by 

k R TA 1 



c+zoo 

dwT{w)z- w 1 (A4) 



c — zoo 

c+zoo \2~w 



~-^(&(3) + A J (21„A + 21„2-l)-iA 3 + 

The first term is the high temperature leading term of the massless limit. Similarly, we find that the asymptotic 
behavior of A\ is 

*bTA i , ,„, ,,2,3 



The most interesting case is the low temperature small mass asymptotic behavior. Define 

ZH = £'(4Ta^- 



Using (|A4[) , we obtain from (|A2[) that 



A 



k B TA 1 



c+io ° Z (w - 2) 

dwT(w)2- w (; R {w + l)- 1 



2nd 2 2m 

Making a change of variables w n- w + 2 and using the identity 

r W = Cr^r^ + 1 



w-2 



we find that 



Similarly, we have 



A 



k B TA 1 
87rid 2 2n~i J c - 



c+?oo 



duT 



w + 3 



Ax 



k B TA 1 
47rld 2 2-7ri 



e+ioo 



+ 2)r 



w + 3 



Cr(w + 3)r ( - ) 2( M ). 



We can consider two cases: $ X <C 1 and A < i? < 1. When t9 A, 

oo ,.oo 



r (I) Z( W ) = J^' dttf- 1 exp {—tn 2 d 2 - tX 2 ) 

n=0 ^° 

~- n •/ 



Thus, when ■& < A < 1, 



ra=0 

fir ^ 

7? I 2 



1440d 3 V 7T 



n— 1 



!5 x2 60. 3 
1 ^A 2 + — A 3 



K 



2imX\ 



Ai 



-K 2 hcA 
' 240d 3 



1-^tA 2 



This comes from the zero temperature term. The thermal correction terms are exponentially small. 
When A < 



n=r u i=0 
oo 



x 2j 



(-1) 



Hence, as A <C $ -C 1, 
A) ■ 



1440d 3 



2 V 2 



1 + ^Cnim 3 ~ ~, 
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l+-a ( lni9-lnA-ln(27r) 



-2j-w 



1\ I? 2 



-4l 



7T 2 ftcA 



240d 3 
up to order A 2 terms. 



3tt 4 
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2. TM contribution 



For the TM contribution, expanding the logarithm in ([TB]) , we find that 



=0 i=l J J V^I+a 5 

Up to first order in a, we find that 



(A5) 



A \ 4j , . z 2 + A 2 - z 2 ( A 



+ Ajaz— 



rz + s/z 2 -\ 2 ) z 2 -\ 2 \z + \J z 2 - A 2 

Therefore, when a <C 1, we have 



where So = -Ao 



ES~Eo + Co + (E,+C,)a, 



n=0 j = l 



2k B TA^,^ f°° , 2 z 2 +A 2 -z 2 / A 



The asymptotic behaviors of Bq have been computed in Section IA 11 Here we compute the asymptotic behaviors of 
Co, B\ and C\. 

First we consider the large mass limit A 1. In the high temperature •& 3> 1 region, the dominating terms are the 
terms with zero Matsubara frequencies. Making a change of variables z > z + A, we find that when A, d 1, 



k B TA™l [°° ( A Y 3 



-2jz 



k B TA [°° , . 4 . / A 



dz(z + A) C= = ) e- 2 ^ 



4?rd 2 J \z + X + Vz 2 + 2z\J (A7) 

^Ae"- fie-" 



4lTd 2 

8^d 2 ' 



Similarly, we have 



B l ~0, 

k B TA 2 _ 2A 
Ci — — 75- A e 

In the region iJ<lcA, 

Co = - — ttt > - / du dzz[ e~ 2jz 



— ^7 — f dz z \l z 2 — A 2 | ; 

ftcA ,3 ,1 



-2jz (A8) 



167T2 d? 
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Similarly, we find that 

127rid 3 

Next we consider the small mass A -C 1 limits. In the high temperature $ 3> 1 region, the dominating term comes 
from the terms with zero Matsubara frequencies. Therefore it is obvious that Bo is exponentially small. For Co, 



k B TA ~ i /•<» / A ^ 4J 
fc B TA y> a 2 r°° C i 



-2j'z 



-2Ajz 



fc B TA 1 
47rd 2 27ri 

Making a change of variables 



/>c+?oo °° />oo / -i \ 4j 

/ ^I»A 2 - Y, 2~ w j~ w ~ 1 / z 1 -™ == 



1 + u 





2>' 
we have 

4j 



^ /"OO 

.7 = 1 



.Er 1 / dv(l -v)(l + v) l ~ w v^-^. 

8 j=l J ° 



The largest pole of this function is at w — — 2. Thus, 







r 






-In 2 - 






2 


4_ 





Similarly, 



C ~-^(A 2 S(0)-A 3 S(-1) + ...) 
And 2 



C kBTA X 3 + 



Next we consider the case where A,i9 <C 1. When t) < A C 1, the leading contribution comes from the zero 
temperature term. As in (|A8p . we have 



Y- dzzV^ii ±_) 



J 

Similar to the case where A < 1 C i?, we find that 

hcA J-^ A 3 ^ / 1 



64^ 



icA 4^ A 3 f 1 , . . , ,■ 5 



247T 2 d 3 
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For C\ , one can show in the same way that it is O (A 4 In A) . 
For we find that as < A < 1, 



hcA 
'3^3 

HcA 

: 37T 2 d 3 



°° poo 

]T / dzz 2 ^z 2 - X 2 e~ 23Z 

OO „qq 

A 4 ^/ dzz 2 V / ^Te" 2; 



1 r c+io ° 



3ir 2 d 3 27ri 
fcvl 1 
127rld 3 27ri 

7T 2 ^cA 



ciw;r( W )2- tu A 4 - tu Ci?.(w) / dzz 2 ~ w Vz 2 -l 



720d 3 



C— ZOO 

C+iOO 



c— lOO 

5 



dw r(w)-^y-2- w \ 4 - w t R (w) 



The case A <C 1? -C 1 is most technical. After some painstaking computation, one can show that up to order A 2 , 
the temperature correction terms are still exponentially small for Co and C\ . For B\ , we have 



2k B TA 00 '<» 

h J ^""" Z2_A 



6 i — ^^2^ <L 



ird 2 
2k B TA c 



/ dwT(w)2- w Q R {w)Y,zl dz(z + z 2 l + X 2 ) 1 ^(2 

J c—iao n=l ^ 



k B TA 1 



7rci 2 27ri 

fc B rA 1 
■nd 2 2iti 



k B TA 1 
7T<i 2 2ixi 

TT 2 hcA 



c+ioo 



dwr(u>)2- ,i 'Ci?(w) 



(A9) 



£>n / dz(z + z 2 )-^ 1 
n=l Jo V 



1 — tO, 2/ 2\-l 

+ A 2 (z + z 2 ) r + 



c+ioo 



dwr(w)2- u 'Cfl(to) 



1/7 — 1 
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720d 3 
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